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Abstract 

Background: Phosphorylation-dephosphorylation cycles (PDCs) mediated by l<inases and pliospliatases are 
common in cellular signalling. Kinetic modelling of PDCs has shown that these systems can exhibit a variety of 
input-output (dose-response) behaviors including graded response, ultrasensitivity and bistability. In addition to 
proteins, there are a class of lipids known as phosphoinositides (Pis) that can be phosphorylated. Experimental studies 
have revealed the formation of an antagonistic kinase-phosphatase complex in regulation of phosphorylation of Pis. 
However, the functional significance of this type of complex formation is not clear. 

Results: We first revisit the basic PDC and show that partial asymptotic phosphorylation of substrate limits 
ultrasensitivity. Also, substrate levels are changed one can obtain non-monotonic bell-shaped dose-response curves 
over a narrow range of parameters. Then we extend the PDC to include kinase-phosphatase complex formation. We 
report the possibility of robust bell-shaped dose-response for a specific class of the model with complex formation. 
Also, we show that complex formation can produce ultrasensitivity outside the Goldbeter-Koshland zero-order 
ultrasensitivity regime through a mechanism similar to competitive inhibition between an enzyme and its inhibitor. 

Conclusions: We conclude that the novel PDC module studied here exhibits new dose-response behaviour. In 
particular, we show that the bell-shaped response could result in transient phosphorylation of substrate. We discuss 
the relevance of this result in the context of experimental observations on PI regulation in endosomal trafficking. 

Keywords: Cellular signalling. Zero-order ultrasensitivity, Phosphoinositide regulation, Endosomal trafficking 



Background 

Biochemical networks have a modular structure [1]. 
The functional modules have different dynamical and 
input-output properties [2,3]. For example positive 
feedback loops can produce bistability while negative 
feedback loops filter noise. An important biochemi- 
cal module in cellular signalling is a phosphorylation- 
dephosphorylation cycle (PDC). Phosphorylation is a 
common post-translational covalent modification of pro- 
teins and lipids, that is mediated by kinases and needs ATP 
to proceed. However, dephosphorylation is mediated by 
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phosphatases and does not need ATP to proceed. Phos- 
phorylation can affect binding properties, localization and 
activity of proteins and receptors [4] . 

Systems with phosphorylation-dephosphorylation 
cycles can exhibit a variety of input-output or dose- 
response behaviors [5]. The level of phosphorylated 
substrate at steady-state is controlled by the kinase- 
phosphatase balance (KPB), i.e., the ratio of total 
active kinase to active phosphatase concentration. If 
the enzymes are far from saturation, the phosphory- 
lated level of substrate is a graded function of the KPB. 
However, if the enzymes are saturated and are oper- 
ate in zero-order regime, an ultrasensitive switch-like 
response is achieved where a small change in the KPB 
can produce a large change in the level of phosphorylated 
substrate [6] . Small modifications to the structure of these 
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biochemical modules, such as the introduction of cooper- 
ativity or product inhibition can significantly affect their 
dynamical properties [7,8]. Cascades of phosphorylation- 
dephosphorylation cycles such as MAP kinase cascades 
also produce ultrasensitivity and amplification [5] . Multi- 
site protein phosphorylation with distributive mechanism 
can produce robust ultrasensitivity [9-11]. Systems with 
multiple phosphorylation sites can exhibit additional 
dynamical properties including multistability [12,13]. 
Multisite phosphorylation can also result in robust 
ultrasensitivity outside zero-order regime through local 
saturation [14,15]. In addition, the order and the distribu- 
tivity of phosphorylation affects the response properties 
of the system [4]. 

During endosome trafficking tagging molecules that 
specify the identity and the fate of different vesicles need 
to be tightly regulated in time. Two classes of molecules 
that have a fundamental role as molecular tags in mem- 
brane traffic are small GTPases and phosphopinositides 
(Pis) [16]. Pis are a type of cellular phospholipids that 
can undergo phosphorylation-dephosphorylation cycles. 
Although Pis are present at relatively low abundance, 
they have important roles in cellular regulation and 
membrane trafficking [17,18]. In particular. Pis can get 
phosphorylated at positions 3, 4 and 5 of their inos- 
itol headgroup, through organelle specific kinases and 
phosphatases [17,18]. Interestingly, several antagonistic 
kinase and phosphatase pairs regulating phosphorylation 
in Pis form complexes (Figure 1). Since these kinsase- 
phosphatase complex pairs have opposing effects, the 
formation of this kind of complexes between them is puz- 
zling and their functional significance is not clear. For 
example, the Vps34 kinase generates PI(3)P on endosomes 
by phosphorylating the third position of the PI inosi- 
tol headgroup and is regulated by VpslS protein kinase 
[19]. Cao et al [20,21] demonstrated that Vps34-Vpsl5 
subcomplex could also interact in vivo and in vitro with 
the direct antagonist MTMl and MTMR2 phosphatases, 
where VpslS mediates the interaction between the Vps34 
and the myo tubular in isoforms. 

Another example of an antagonistic kinase-phosphatase 
pair in yeast and mammals is the Fabl PI(5)P 5-kinase 



(also called PIP5K3 or 'PIKfyve') and the Fig4 PI(3,5)P2 
5-phosphatase (also called Sac3) [23-25]. Moreover, the 
existence of two new kinase-phosphatase complexes 
which play a role in the regulation of PI(5)P has also been 
hypothesized [22]. In particular, it has been proposed 
that a PIK3C (class-I PI 3-kinase) and myotubularin com- 
plex regulates the interconversion between PI(3,5)P2 and 
PI(5)P [22]. Less is known about the PIP5K2 4-kinases and 
PI(4,5)P2 4-phosphatases which regulate the interconver- 
sion between PI(4,5)P2 and PI(5)P, however, the existence 
of a corresponding kinase-phosphatase complex has also 
been speculated [22]. The identified and hypothesized 
kinase-phosphatase pairs are summarized in Figure 1. In 
summary, there is a growing evidence that a number of 
kinases and phosphatases can form a complex, but the 
exact role of such a kinase-phosphatase duo remains to be 
investigated [26]. Moreover, there is some evidence that 
this kind of complex formation is evolutionarily conserved 
which suggests functional significance [22]. Appropriate 
mathematical modelling may shed light on the functional 
roles of these complexes. 

In this paper, we use mathematical modelling to inves- 
tigate the dose-response in an extended phosphorylation- 
dephosphorylation cycle that includes complex formation 
between the kinase and phosphatase. In particular, we 
ask under what conditions the PDCs can exhibit a non- 
monotonic dose-response. To this end, we first revisit the 
basic PDC (Figure 2a) and show that depending on the 
parameters of the system it is possible to only achieve par- 
tial phosphorylation of the substrate even for large KPB. 
We then use this results to refine estimates of zero-order 
ultrasensitivity. We further show that at a fix kinase- 
phosphatase balance, if we change the substrate levels, 
it is possible to obtain non-monotonic bell-shaped dose- 
response. However, this behavior is only obtained over a 
narrow range of parameters. 

To investigate the properties of PDC with complex for- 
mation, we need to assume enzymatic activity of the 
kinase-phosphatase complex. Thus, we study four dif- 
ferent variants of the extended PDC module, where 
the complex exhibits kinase and/or phosphatase activity 
(Figure 2b-e). We observe that the module can produce an 
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Figure 1 Kinase-phosphatase complexes. Identified or liypotliesized l<inase-pliospliatase complexes in the PI regulation system. Full line 
represents in vivo and dashed line in vitro evidence. The data is based on the reference [22]. 
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Figure 2 Reactions diagrams, (a) Reaction diagram for tine basic 
model of pliospliorylation-depliospliorylation cycle, (b-e) Reaction 
diagrams for the extended model with kinase-phosphatase complex 
formation. The complex can have kinase and/or phosphatase activity. 
The grey letters in the enzymatic complex name correspond to the 
inactive enzymes. There are four variants of the extended model with 
enzymatic complex having both inactive K and P (b), active P and 
inactive K (c), active K and inactive P (d) and both active K and P (e). 



ultrasensitive response outside the zero-order regime. In 
addition, we observe robust bell-shaped non-monotonic 
dose-responses in a variant of the module where the 
complex has phosphatase activity. Finally, we discuss our 
results in the context of PI regulation and endosome 
trafficking. 

Results and discussion 

Basic model of phosphorylation-dephosphorylation cycle 

We shall consider a system in which a substrate can exist 
in the unphosphorylated form S and phosphorylated form 
Sp, The substrate could represent a protein or a PI. A 
kinase K and a phosphatase P mediates phosphorylation 



and dephosphorylation (Figure 2a). S and Sp are controlled 
by the following enzymatic reactions: 



Sp+P 



^2 



SpP 



S^P 



■ K ^ SK Sp + K, 

^-2 



(1) 



(2) 



where, the rate constants are denoted A.i,A.2 for binding 
an enzyme and substrate, A._i, A._2 for dissociation of the 
complex and /q, k2 for the catalytic reaction. The associa- 
tion rates (A./) have units /xM~^s~^, whereas dissociation 
(A._/) and catalytic (Iq) rates have units s~^. The kinetic 
equations governing the time evolution are 

[^p] = ^-1 M + f<2[SK] -Ai [Sp] [P] (3) 



dt 

^ [SpP] = + ki) [SpP] + ki [Sp] [P] 



dt 
d 

dt 



[S] = k^2[SK] +ki [SpP] -k2[S] [K] 



[ SK] = - (A_2 + /:2) [ SK] +k2 [ S] [ K] . 



(4) 



(5) 



(6) 



This system is complemented by conservation equations 

Kt=[K]+[SK]. Pt=[P] + [SpP]. (7) 

St=[S] + [Sp]+[SK] + [SpP], (8) 

where the total enzyme /<Q, Pt and substrate St concentra- 
tions are in /xM and are often assumed to be constants. At 
equilibrium. 

Ml [SpP] = [Sp] [P] , M2[SK] =[S] [K] , (9) 

where the Michaelis constants Mi = (i = 1, 2) are 

in /xM. 

The system (3)-(8) has been extensively studied over 
the past 30 years since the key study of zero-order ultra- 
sensitivity by Goldbeter and Koshland [6]. Most studies 
look at the steady-state dose-response of the average 
phosphorylated substrate versus the kinase-phosphatase 
balance (Kt/Pt)^ As a measure of phosphorylated sub- 
strate one can either look at the fraction of total phos- 
phorylated substrate {R = ([Sp] + [SpP])/St) or at the 
fraction of free phosphorylated substrate {R^ = [Sp] /St), 
The relevant functional form could be either free or total 
phosphorylated substrate depending on the specificity 
of substrate binding sites for phosphatases and recog- 
nized targets. In this study, we follow other works that 
use the fraction of total phosphorylated substrate [R) to 
study the dose response in the PDCs [7,27]. An important 
approximation used to analyse these systems is the quasi 
steady-state assumption (QSSA) [28] which assumes that 
the enzyme-substrate complex concentration remains 
approximately constant over time (apart from a fast initial 
transient). Under QSSA, total {R) and free {R') follow the 
same trends. 



Szomolay and Shahrezaei BMC Systems Biology 201 2, 6:26 
http://www.bionnedcentral.conn/l 752-0509/6/26 



Page 4 of 14 



Most studies assume that the kinase-phosphatase 
balance is changed by varying the level of kinase {Kt)> 
while keeping the phosphatase levels (Pt) constant. This 
is probably the most relevant regulation that happens in 
cellular systems. For example, in a MAPK cascade the 
activity of upstream kinase is increased while effectively 
the level of phosphatase remains constant. However, this 
way of varying KPB is asymmetric with respect to R, This 
is why our analysis shows that for low KPB the fraction 
R is always close to zero (no phosphorylation of sub- 
strate), however for large KPB the fraction R can take 
any value smaller than one (partial asymptotic phospho- 
rylation of substrate). There are two other ways that KPB 
can be changed. Another asymmetric method is to keep 
kinase levels constant and change the level of phosphatase. 
An example of this kind of regulation has been recently 
reported in the regulation of yeast mating pathway [14]. 
A symmetric way of changing KPB is to keep the total 
level of kinase and phosphatase constant and change the 
ratio only. This is probably less relevant to cellular sys- 
tems. However, this way of changing KPB will always 
result in full asymptotic phosphorylation of substrate at 
large KPB. Finally, it should be noted that while over some 



signalling time scales KPB might increase monotonically, 
over longer time scales there are typically negative feed- 
backs in place that would decrease KPB, e.g., by increasing 
the levels of phosphatases [29,30]. In this study, we focus 
on short-term monotonic changes of KPB ignoring the 
slow negative-feedback processes. 

Partial asymptotic phosphorylation of substrate and Hill 
numbers in the basic model 

Under QSSA, the asymptotic limit of R for large kinase- 
phosphatase balance (e.g., as Kt oo), denoted by R^l 
is always one, meaning that the substrate gets fully phos- 
phorylated. However, this is not the case in general when 
QSSA is not valid, as shown using the exact steady-state 
solution of (3)-(8) (see Methods). In fact, the asymptotic 
phosphorylation can be very small for some parameter 
values as seen using the analytically calculated R^^ depen- 
dence on the relevant parameter ratios St/Pt>o(, ki/k2 
in Figure 3(a). Intuitively, the substrate is only partially 
asymptotically phosphorylated when the phosphatase can 
sequester a significant fraction of substrate by either satu- 
rating the substrate (small ^), having higher catalytic rate 

(small a = ^), or having a higher affinity for substrate 
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Figure 3 and nn for the basic model, (a) is computed analytically as St/Pt, a and X] /Xj is varied. Baseline parameter values used for the 
basic model, are: = ] s-\ X^ = 1 M"^ s-^ /c] =1 s"^ , A._2 = 1 s"^ , X2 = ^ /xM"^ s"^ , /c2 = 0.01 s"^ , Pt = ^ /xM, St = ^0 /xM. (b) nn is 
computed analytically as St/Pt, a and X] /X2 is varied. Baseline parameter values used for the basic model are: = ] s~\ X] =1 /xM~^ , k] = 
1 s-^ A._2 = 1 s-^ X2 = 1 /xM-^s-^ k2 = l s-^ Pt = l /^M, St = 100/xM. 
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than kinase (large ^) (Figure 3(a)) [31]. Similar results 
have been previously obtained for other protein modi- 
fication cycles [32]. The observed limits in Figure 3(a) 
are r'^ ^ 1 for St/Pt CO and 1^ a/(a + 1) for 
^i/^2 oo; these results are obtained in the Methods 
section. 

Goldbeter and Koshland have introduced analytical for- 
mulas for the response coefficient (which is defined as 
the ratio of the kinase concentration required to give 
90% response relative to the concentration required to 
give 10% response; assuming constant phosphatase con- 
centration), both using QSSA and the exact steady-state 
solution. These results, however, do not take into account 
the limit of R which in general can be smaller than one 
and based on which the 10% and 90% value should be cal- 
culated. Here, we have used results for R obtained analyt- 
ically to introduce a new response coefficient formula and 
hence, to find the Hill-number (see methods). To illus- 
trate our analysis, we have looked at the dependence of the 
analytically calculated hh with respect to relevant param- 
eter ratios St/PtyOiyXi/)^2 (Figure 3(b)). Comparison with 
numerically obtained values of R^^ and hh (obtained by 
fitting a Hill function; not shown) confirms our analyti- 
cal results. As expected, hh is very large at high substrate 
to phosphatase concentration ratio (zero-order regime; 
Figure 3(b)). Interestingly, large changes in hh are also 
obtained when a and A1/A2 are varied. Comparing results 
in Figure 3(a) and Figure 3(b) suggests that high Hill 
numbers are only achieved when is close to zero. 
However, modest ultrasensitivity is present even at very 
low asymptotic phosphorylation levels of substrate. These 
results are also consistent with another work which sug- 
gests that substrate sequestration can significantly reduce 
ultrasensitivity [31]. 

Bell-shaped dose-response in the basic model 

Most biochemical signalling that involves phosphoryla- 
tion is mediated by changes of KPB by either regulat- 
ing the level of active kinase or active phosphatase as 
explained above. However, one can imagine cases where 
the substrate level can change, while the KPB is relatively 
constant. Therefore, it is interesting to ask what kind of 
dose-response one can expect when the substrate levels 
are changed. 

To address this question we introduce asymptotic phos- 
phorylation of substrate for large substrate levels (as Kt 
00) denoted by 7?oo« As shown in the Methods, when 
the maximum speed of kinase activity is greater than the 
maximum speed of phosphatase activity (^ < Kt), the 
asymptotic phosphorylation is i^oo = 1> since kinase activ- 
ity wins at large St. For this case we observe a monotonic 
increasing dose-response with mild ultrasensitivity with 
Hill coefficients about 2 (Figure 4(a)). More interesting 
behaviour is achieved for the case when the maximum 



speed of kinase activity is less than the maximum speed of 
phosphatase activity (^ > Kt) and the asymptotic phos- 
phorylation is i?oo = 0- In this case, there are specific 
parameters, where a bell-shaped non-monotonic dose- 
response can be observed (Figure 4(a)). As St is increased 
from zero, R is initially a constant that is set by the KPB. 
For large St, it eventually approaches its asymptote of zero. 
For most parameters, the transition from constant R to 
zero is monotonic (Figure 4(a)). However, as shown in 
Figure 4(b), over a narrow region of parameters above the 
line Pt /Kt = a, it can exhibit a significant non-monotonic 

bell-shaped response. The bell-shaped response is pro- 

s s 
duced as the transition from R^ = 0 to R^ = 1 is 

approached at Pt/Kt = a (See Eq. 26 in the Methods). 

Since, this behaviour is only achieved with fine tuning 

of the parameters and therefore is not expected to be of 

biological relevance. 

Model with kinase-phosphatase complex formation 

In this section we consider an extension to the basic model 
of phosphorylation-dephosphorylation cycle (Figure 2). 
We assume that the phosphatase and kinase can form 
a single protein complex PK given through a reversible 
reaction (10): 

P-\-K ^ PK (10) 

K-l 

Depending on which enzyme is 'on' or 'off within the 
kinase-phosphatase complex, the substrates S/Sp could 
be further phosphorylated/dephosphorylated by PK as 
shown in (11) and (12): 

Sp+PK ^ SpPK X S + PK (11) 

S-\-PK ^ SpPK % Sp+PK (12) 

In (10)-(12), the parameters a:i,A3,A4 and a:_i,A_3,A_4 
correspond to the association/dissociation rates, respec- 
tively, and k'^y represent the catalytic reactions. Here we 
define a parameter co = k-\/ki, similar to the concept 
of a 'dissociation constant', which expresses the affinity of 
binding between the phosphatase P and kinase K. 

There are four possibilities as shown in Figure 2; both 
enzymes in the complex could be enzymatically inac- 
tive (Figure 2(b)), only phosphatase active (Figure 2(c)), 
only kinase active (Figure 2(d)) or both kinase and phos- 
phatase active (Figure 2(e)). This system is described by 
the following kinetic equations 

^ [Sp] = k-i [SpP] + k2[SK] -ki [Sp] [P] 

-\-X-3[SpPK]-X3[Sp][PK]+kdSPK] (13) 
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Figure 4 Bell-shaped dose-response for the basic model, (a) Bell-shaped dose response for the basic model when St is varied. Parameters used 
are: X-^ = X^ = = X-2 =0.]s-\k^ = ] s-\ k2 = 0.] s"^ , Kt = 2000 /xM, En = 0.1 001 (-), En = O.IE^ ( — ), En = 0.08Ek (-.) (b) 2D 



parameter plot when a and Pt/Kt are varied. Parameters used are A._ 



:A2 



: 0.1 s-\a = 0.1, Kt = 2000 /iM. R^rliax denotes the 



maximum of /? with respect to St and Rq = limSf^o+ /^.Within a narrow region above the line Pt/Kt = a the ratio /^max/^o > 1-05, i.e., there is at least 
a 5% increase in /?max/^o compared to Rq. 



^ rc Dl n ^/ ^rc Dl^. rc iFDi ^A^ -[SK]= + k2)[SI<]+k2[S][I<] (16) 

— [SpP\ = -{X-i + ki) [SpP\ + Xi [Sp\ [P] (14) dt 

^ ^ [ PK] =Ki[P][ K] -K-i [ PK] +(A_3 + ^3) [SpPK] 
-[S]= X.2ISK] +k, [SpP] - X2[S] [K] +ks [SpPK] _ ^5^] [^^j ^ ^^^[^^^j _ 

+ X-^{SPl<\-X^{S\{PK\ (15) -X^{S\{PK\ (17) 
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- [SpPK] = As [Sp] [PK] -(A_3 + ks) [SpPK] (18) 



- [ SPK] = X4S][ PK] - (A_4 + h) [ SPK] (19) 
at 

coupled with conservation equations 
[P] + [SpP] + [SpPK] ^[PK] = Pu 

(20) 

[K] +[SK] +[SPK] +[PK] = Ku 

[Sp] +[ ^] + [SpP] +[SK] + [SpPK] +[SPK] = St, (21) 

We note that the the extended model (13) -(21) is too com- 
plicated to find the corresponding steady-states analyti- 
cally. In the following we have used numerical integration 
of the kinetic equations above and variants of the QSSA 
to study the dose-response in this system. The ration of 
phosphorylated substrate in the extended model is now 
defined 2iS R = ([Sp] + [SpP] + [SpPK])/St, Our goal is 
to determine the dose-response for as a function of KPB 
in the extended model. We have mainly focused on the 
kinase active case and the phosphatase active case which 
are the two extreme scenarios (Figure 2(c)-(d)). 

Ultrasensitivity 

We observe ultrasensitivity for both the kinase active- 
phosphatase inactive and the phosphatase active-kinase 
inactive cases outside the zero-order regime, where St < 



Pt (Figure 5). For the phosphatase active case, the mecha- 
nism is similar to the previously described reversible com- 
petitive inhibition, where the substrate and the inhibitor 
cannot bind the enzyme at the same time [5]. In this case, 
when the level of K is increased, initially it binds to P and 
forms KP that acts as a phosphatase and R remains close 
to zero. However, at some level of /<, most of the free P 
will be in complex with IC Above this level a small increase 
in K changes the phosphorylation level of substrate {R) 
nonlinearly, producing ultrasensitivity. 

The mechanism of ultrasensitivity for the kinase active 
case is novel. In this case, when the level of K is increased, 
it binds to P and forms KP that acts as a kinase. Therefore, 
increasing levels of K has the compound effect of increas- 
ing the kinase activity and at the same time reducing the 
phosphatase activity, resulting in a nonlinear increase in R 
and an ultrasensitive response. As shown in the Figure 5, 
both ultrasensitivity mechanisms are independent of the 
enzyme saturation and can produce sharp responses in the 
parameter regime where the basic model cannot. 

Bell-shaped dose-response curves 

The active phosphatase case, when the kinase is off and 
the phosphatase is on in the protein complex is inter- 
esting since it can exhibit a non-monotonic response. 
Suppose that the affinity between P and K is low and 
the PK complex is a stronger phosphatase than P alone 
(with a lower i?^). In this case, when the level of K 
is increased initially, the KP production is small. How- 
ever, the kinase activity soon takes over the phosphatase 
activity and R is increased close to 7?^. As K is increased. 




Figure 5 Ultrasensitivity. (a) R, total K, P and KP are numerically calculated as KPB is changed in the extended model with only kinase active 
(Figure 2d). Parameters used are: = 1 0 s~\ A.] = 1 /xM~^ , k] = 1 00 s~\ A._2 = 1 0 s~\ A.2 = 1 /xM~^ s~\ /c2 = 1 00 s~\ A._4 = 1 0 s~\ 
A.4 = 100/xM-^s-^ /C4 = 0.01 s-\ K-] =]Os-\k] =] /xM-^s-^ Pt = 10 /xM, St = 1 /xM. Hill coefficient for basic model and active /( is n^^ = 1 
and hh = 1 .8, respectively, (b) R, total K, P and KP are numerically calculated as KPB is changed in the extended model with only phosphatase active 
(Figure 2c). Parameters used are: =10 , h = ] /xM~^ s~\ /c] = 1 00 s~\ A._2 = 1 0 s~\ A.2 = 1 /xM~^ , /c2 = 1 00 s~\ A._3 = 1 0 s~\ 
A.3 = 1 00 /xM-^ s-\ k3 = 0.01 s-^ , K-] =]0s-\ K] =] /xM"^ s-\ Pt = ^0 /xM, St = 1 /xM. Hill coefficient for basic model and active K is 
riH = 1 and nn = 1 .95, respectively. In both case the bottom panels show that the level of enzymes remains above the substrate, so we are far from 
the zero-order regime. 
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Figure 6 Bell-shaped dose-response in the active phosphatase 
case (Figure 2c). (a) Bell-shaped response in the zero-order regime 
when St > Pt. Parameters used are: = ]0s~\ X] =1 /xM~^ 

/Cl = 1 S-\ X-2 = 10s-^ X2 = 20/xM-^S-^ /C2 = 1 s-\X-3 = 

]0s-\ A.3 = 1 fzM-^s-\ /C3 = 100cs-^ = 50s-\ /C] = 
1 ixM-h-\ Pt = 1 /xM, St = 10/xM. The Hill coefficient for the 
basic model is nn = 2.1 5. (b) Bell-shaped response outside the 
zero-order regime when St < Pt' X- 
. X-2 = ^0s-\X2 

-1c-1 



lOs-^ 

-1c-1 



-1 



,X] = 1 fzM 

/(] = 1 s-\ X-2 = ^0s-\X2 = 20/xM-^s-\ /C2 = 1 s-\X-3 = 
10s-^ ^3 = 1 /xM-^s-^ /C3 = 100s-^ = 50s-\K^ = 
1 fzM-^ s-\ Pt = ^ /xM, St = 0.1 fzM. The Hill coefficient for the 
basic model is n^^ = 1 . In (a)-(b), the full line corresponds to the basic 
model, the dashed line is numerical results for the extended model 
and the dotted line shows the tQSSA approximation. 



the KP production becomes more dominant, and, since 
KP is a stronger phosphatase, the is dropped to a 
lovv^er level, hence, producing a bell-shaped response. As 
Kt oo, the limit of the response curve has been derived 
earlier for the basic model. This measure can give an esti- 
mate of the peak level of R, However, this -formula 
could be used for the active P case as well. Indeed, at 
high Kt concentrations the free phosphatase is close to 
zero and we effectively reduce to the basic model with 
only K and PK around. Hence, replacing the parameters 
A.1, A._i, /q by A.3, A._3, /cs in (27) gives the desired result for 
the asymptotic level of R, 

Bell-shaped responses could be obtained both at high 
enzyme saturation {St ^ Pt)> and at low enzyme satu- 
ration (St ^ Pt)' Figures 6(a)-(b) compare the response 
curve for the active P with the response curve for the 
corresponding basic model for both scenarios. The only 
difference is that for the low enzyme saturation we have 
an ultrasensitive rise to the peak of the bell-shaped 
response, whereas for the high enzyme saturation there is 
a linear rise to the peak. 

The steady-state solutions of the extended model can 
only be defined implicitly. Hence, we can use the total 
quasi-steady-state-approximation (tQSSA), which is 
more generally valid than the QSSA [33,34]. The tQSSA 
method has been successfully applied to a model with a 
pair of enzymes and substrates (analogous to the basic 
PDC) [27] and to coupled PDCs [35]. The details of the 
tQSSA method are explained in the Methods section. The 
results of tQSSA calculation is consistent with the numer- 
ical results in Figures 6(a)-(b). The fact that tQSSA but not 
QSSA can explain the bell-shaped response suggests that 
the dynamic of enzyme-substrate complex contributes to 
the presence of the non-monotonic response. 

We have investigated the robustness of the bell-shaped 
response by varying some of the key parameters. We 
have characterised the bell-shape response by looking 
at the ratio of the maximum response to asymptotic 
phosphorylation as a function of dissociation constant of 
the phosphatase-kinase complex a; = the Michaelis- 

Menten constant M3 = ^"^^^"^^ of the complex and the 
total phosphatase level Pt (Figure 7). It is possible to 
obtain significant bell-shaped responses over a large 
range of parameters. The response has the highest relative 
peak for large values of co and Pt (Figure 7(a)), however 
for each parameter set there is an optimal value of M3 
(Figure 7(b)). Also, we have compared the transient time 
needed to approach steady-state in the extended model 
with active P (Figure 8). We observe that the transient 
time for the extended model follows closely the transient 
time for the basic model and it goes through a maxi- 
mum at the threshold of ultrasensitivity. This is because 
the enzymatic complex formation is fast compared to 
phosphorylation and dephosphorylation. 



Possible role for complex formation in phosphoinositides 
regulation 

We propose that complex formation can contribute to 
the transient regulation of Pis observed experimentally 
in endosome trafficking [36]. Suppose the kinase is being 
recruited to the vesicle with a constant rate (slower than 
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Figure 7 Robustness of the bell-shaped response. To illustrate the robustness of the bell-shaped response, the ratio of R^ax> the maximum of R 
with respect to Kt, to is calculated as multiple parameters are varied. The red regions correspond to parameter values with the largest 
bell-shaped response, (a) co and Pt are varied as a variable in co). (b) co and M3 are varied as a variable in co and A.3 as a variable in M3). 
Parameter values used in the extended model with only phosphatase active in (a) and (b): A._i = 1 0 s~^ A.] = 1 /xM~^ s~\ k] = ] s~\ k-2 = 
]0s-\ Xj = 20/xM-^s-\ /C2 = 1 s-\X-3 = 10s-\ X3 = 1 f^M'^ s-\ /C3 = 100s-\ = 10s-\ K^ = 1 fiM'^ s-\ Pt = 1 /xM, St = 10/xM. 



the time scale for the transient behaviour in Figure 8) and formation and the complex has phosphatase activity, in 
there is a fixed level of phosphatase on the vesicle. This the parameter regime where we have bell-shaped dose- 
produces a gradual change in KPB. If there is complex response as discussed above, the level of phosphorylated 
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200 




Figure 8 Transient times to approacii steady-state for the basic 
and extended model. Parameter values for the basic model: = 
]0s-\ X^ = 1 /xM-^s-\ /ci = 1 s-\X-2 = 10s-\ = 20 /xM"^ 
s~\ kj = 1 s~\ Pt = 1 /xM, St = 10 /xM. Parameter values for the 
extended model: = ]0s~\ X] = ] /xM~^s~^ k] = 1 s~^A._2 = 
]0s-\ Xj = 20/xM-^s-\ /C2 = U-\ K-^ = 50s-\ K^ = 1 /xM"^ 
s-^A._3 = 10s-^ A.3 = 1 /xM-^s-^ /cs = ]00s-\Pt = 1 /xM, St = 
10/xM. The transient times for /?= ([Sp] + [SpP] + [SpP/(])/Sr were 
estimated such that the absolute difference between the positive 
steady-state and the time solution at the transient time is within a 
small tolerance and that the transient times do not change too much 
if the tolerance is varied. Hence, we chose the tolerance to be 1 0~^. 



PI goes through a maximum as KPB is changed over time. 
This can produce a fast and accurate transient change in 
the phosphorylation of PL 

Botelho [37] has speculated on the possible role of enzy- 
matic complex formation in transient regulation of Pis but 
he has proposed an alternative mechanism to achieve a 
transient behaviour. This alternative explanation assumes 
that there is a signaling trigger that activates the phos- 
phatase relative to the kinase in the enzymatic complex 
leading to the down-regulation of phosphorylated PI on 
the membrane. Therefore enzymatic activity of the com- 
plex switches from kinase to phosphatase upon signalling. 
The signalling trigger might be the recruitment of other 
vesicles or GTP hydrolysis by GTPases [37]. The Fabl- 
Figure 4 and Vps34-MTM complexes in PI regulation 
could follow this concept (Figure 1). 

In Figure 9 we have compared the time evolution of 
phosphorylated substrate (PI) following the two alterna- 
tive mechanisms mentioned above. In the first scenario, 
we assume a hypothetical recruitment of kinases over 
time to the membrane, fixed levels of phosphatase on 
the membrane and the active phosphatase activity in the 
PK complex. In the second scenario following Botelho 
[37], we assume a fixed level of kinase and phosphatase 
and we assume a switch in the enzymatic activity of the 
PK complex from kinase to phosphatase at ^=300s. Both 
scenarios can result in fast transient phosphorylated sub- 
strate which peaks at 300 seconds, while our proposed 



mechanism results in faster decay as Kt increases over 
time. While both phosphorylation profiles in Figure 9 are 
in qualitative agreement with the observed profiles in the 
process of endocytosis [36], however, the dynamic of avail- 
able kinase is quite different in the two models and can 
be used to experimentally test our proposed hypothesis. 
Note that in our proposed mechanism, the asymptotic 
phosphorylation level of the substrate does not get to 
zero, however, by increasing the phosphatase activity of 
the complex, one can make this asymptotic phosphoryla- 
tion arbitrarily small. Also, a counter intuitive aspect of 
our proposed model is that if the kinase levels are down- 
regulated, then the phosphorylation of the substrate could 
even increase again. So, one requires no down-regulation 
of the upstream kinase over the relevant time-scales of the 
transient PI profiles. Ultimately, the response of the PDC 
module should be considered in the context of the larger 
PI regulation network (Figure 1) to be able to obtain a full 
understanding of the role of enzymatic complex formation 
in the endosome trafficking. 

Conclusions 

In this paper, we have studied the steady-state dose- 
response in basic PDC motif and a modified PDC motif 
using analytical and numerical methods. We first showed 
that there are parameter regimes where even at large KPB 




time (seconds) 

Figure 9 Transient phosphorylation of phosphoinositides due to 
enzymatic complex formation. Comparison of the two scenarios; 
solid line represent the case when a signalling trigger is simulated. 
Parameters used are the following: A._i = ]0s~\ X] = ] /xM~^ , 
/ci = 1 s-\ A._2 = 10s-\ X2 = 1 ixM-h-\ kj = 1 s-\/c_i = 
50s-^ /ci = 1 Kt = 0.1 /xM, Pt = 0.2 /xM, St = 300 /xM. 

The other parameters are A._4 = 1 s~^ A.4 = 75 /xM~^ s~\k4 = 
lOOs-^ fort< 300sandA._3 = 10s-^A.3 = 10/xM-^s-^/c3 = 
500 , for t > 300 s. Dashed line represent the case when kinases 
are recruited linearly in time as Kt = 0.0003t. Other parameters used 
are the following: = ]0s~\ X] = ] ixM~^s~\ k] = ] s~\ 
X-2 = 10s-^ X2 = 1 /xM-^s-^ /C2 = 1 s-\ = 50s-^ /ci = 
1 /xM-^s-^ = 10s-^ A.3 = 10/xM-^s-^ ks = 500s-^ Pt = 
0.2 /xM, St = 300 /xM. 
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the phosphorylated ratio of substrate remains small. This 
is not expected from the QSSA analysis of the basic PDC, 
and it is a consequence of efficient phosphatase binding 
to a substrate. This results in modified estimates of the 
Hill number in zero order ultrasensitivity. In addition, we 
observe that partial asymptotic phosphorylation of sub- 
strate limits the ultrasensitive behavior. Over a narrow 
parameter regime, one can obtain a bell-shaped dose- 
response when the substrate level is changed for a fix KPB 
level. We then studied a modified PDC motif where the 
antagonistic kinase and phosphatase can form a complex. 
Depending on the activity of the kinase and phosphatase 
within the complex, the modified system has four dif- 
ferent cases. We investigated the possible dose-response 
curves in all scenarios. We observed a novel form of 
ultrasensitivity arising from complex-formation that can 
function outside the zero-order regime. Also, we report a 
non-monotonic bell shaped dose-response for the phos- 
phatase active-kinase inactive case that can function over 
a wide-range of parameters. 

Complex formation between antagonistic enzymes have 
been reported in the PI regulatory network (Figure 1). 
Motivated by the observed bell shaped dose-response, 
we propose the complex formation between the kinase 
and phosphatase may contribute to the observed tran- 
sient regulation of Pis during endosome maturation [36] . 
It has been recently proposed that the complex forma- 
tion observed in several kinase-phosphatase pairs in the 
PI regulation may enhance the temporal regulation of 
the Pis [37]. Here we have shown that if the complex 
has enhanced phosphatase activity and negligible kinase 
activity, then the dose-response will be bell-shaped. This 
mechanism produces a transient increase in the phospho- 
rylation level of the Pis as enzymes are recruited to the 
endosomes and PKB is changed over time. There is no 
need for additional regulators or triggering mechanisms, 
but only monotonic change of the PKB on the endosome 
by enzymatic recruitment. Experimental measurement of 
enzymatic activity within the complexes are still miss- 
ing. Our model predicts that an enhanced phosphatase 
activity within the complex gives rise to bell-shaped dose- 
response. Also, we predict that alteration of the total level 
of enzymes or the affinity of enzymes in the complex will 
have a significant effect on the height of the bell-shaped 
response. 

Our study illustrates that simple alterations of the reg- 
ulatory motifs in biochemical networks can have a signif- 
icant consequence on the response and functionality of 
these systems. Signalling systems are typically composed 
of several interlinked regulatory motifs. To ultimately 
understand the function of the whole system one needs to 
put the motifs together and investigates the coordinated 
dynamic behaviour of all the components as they interact 
with one another. In the context of the PI regulation, as 



shown in Figure 1, there are several PDC motifs with enzy- 
matic complex formation. As shown here, these motifs can 
produce ultrasensitive or non-monotonic dose-response 
depending on their biochemical parameters. The chal- 
lenge that remains is to quantify these parameters and 
explain the overall function by investigating the interplay 
between the dynamic of individual motifs. 

Methods 

Limit of the average fraction of phosphorylated substrate 

for Kt^oo and Sf ^ oo 

We derive limiting values for the average fraction of phos- 
phorylated substrate R = ([Sp] + [^pP])/^^ at steady-state 
as Kt oo and St oo, respectively. First, we will 
implicitly define the steady-state [SK] as a solution of the 
cubic equation 



[SKf +Qi[SKV +Q2[SK] +Q3 = 0. 



(22) 



The coefficients Q/ (/ = 1,2,3) are of the form Qi = 
-Kt + vi, Q2 = uiKt + V2, Q3 = U2Kt with 

PtSt 



Pt , Mi+St 

ui = — + -—, 

a Of + 1 



U2 = - 



Qf(Qf + 1) 



Pt Mi+M2 + St Pt(Kt + St) 

Vl = — , V2 = — — 

a Of + 1 a{a + 1) 

where a = ^ . It follows from (3) and (9) that a satisfies 

a[SK]=[SpP\ 
and hence, 

Pt 



(23) 



[SK]< 



(24) 



for Pt > 0. The asymptotic expansion of [ SK] with respect 
to powers of ^ in (22) implies that the asymptotic limit 

of [SK]f i.e., 5/<oJ = limKt->oo[SK] satisfies the quadratic 
equation 

(SI& - UlSK^ - M2 = 0. 

By condition (24) we obtain that 



Ml - yjul + 4M2 



(25) 



Note that Ml = f + and M2 = - Jfiy do not 

depend on M2 and hence, SK^ does not change when X2 
or X-2 is varied. Intuitively, this is because for high kinase 
levels the substrate will be saturated and the SKo<^ will 
not be influenced by the dissociation constant of kinase- 
substrate binding. By (9) and (23) the phosphorylated 
substrate fraction R could be re-written as 



St \Pt-a[SK] ) 



(26) 
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For Kt oo, the limiting value of R denoted by 
becomes 



where is given by (25). It is shown in the Methods 
section (Corollary 1) that 7? < 1 as expected. 

For St oo, we also include a statement about the 
limiting value of R denoted by R^, which will be useful 
later on (for proof see Corollary 2 in the Methods section). 
If f < Kt, then [SK] ^ ^ and ^ 1 as 5^ ^ oo, 
otherwise, [ SK] Kt and ^ 0 as 5^ ^ oo. That is. 



Ml 



■ aSK 



■Kt 



+ 1 



(27) 



r: 



St 



1 if ^ < Kt 
0 if ^ > Kt. 



Corollary 1-3 
Corollary 1. The ratio R satisfies 



R. 



St 



+ 1 < 1. 



Proof It is sufficient to show that < 1 since ^ = J^h = 

d^]^^^ > 0 and hence, R is monotone increasing 

with respect to Kt and R^^ is its maximum. R^^ < 1 is 
equivalent to the condition 

a^iSK^f - a(Mi + Pt + St)SK^^ + StPt > 0. (28) 



To show that (28) holds we use that (SK^^f - uiSK^ - 
U2 = 0 and that Pt > aSK^ (this follows from (24)). 
Multiply (28) by^^to get 



{SK'^f-[ui- 



\SK^-U2>0. (29) 



Note that the left-hand side of (29) is equivalent to 



SK^t [ ^ 



SK'' 



and this concludes the proof. 



> 0 



^ andR 



□ 

-> 1 
0 as 



Corollary 2. If ^ < Kt, then [SK] ■ 
as St ^ oo.If^ > Kt, then [ SK] Kt and R 
St oo. 

Proof Denote SK^ = limSf^ool SK]. It follows from 
(22) that SK^ satisfies the quadratic equation 

C3(SK^f + (-csKt + ci)SK^ - ciKt = 0 

and then 



sk!^ 



(-Kt-^)±\-Kt^ 



Ell 



Thus, SK^ = ^ if^ < Kt and SK^ = Kt otherwise. 



Obviously, if SK^ = Kt, then R% = 0. To see that if 



SK^ 



oo. Indeed, note that [S] = ^t-iSK] 



1, we show that [ Sp] St as 

M2Pt 



aKt-Pt 



as 



^, then R^c^ 

St 

St oo. Hence, [Sp] = St-[SK] - [SpP] -[S]^ 5^ as 

s p 
St ^ oo and so 7?oo = 1 when -f < Kt. □ 



The following Corollary is introduced without its proof, 
although, the arguments would be similar to Corollary 2. 

Corollary 3. It holds that: 

Pt 



i) SK'^ - 

ii) SK^^ 



as St 
0 asa - 



oo and R^ 
oo and 



< then SK; 



las Xi OQ.If^ 
andl^ -^-r aski 



■Kt 
oo 



> then SK, 



oo. 



1 as St ^ oo 

^ I as a ^ oo 
as Xi ^ oo and I^^ 

■Kt 



OO 



Hill numbers 

The Hill number is defined as 



log(81) 



logiK^/Kl') 



where Kf^/K}^ is the response coefficient [6]. Solving for 
[ SK] from (26) and by condition (24) we obtain 



[SK]: 



Ml +Pt+RSt - ^(Mi+Pt + RSt)^-4.RStPt 



2a 



We also find Kt from (22) to be 

[SK]^^vi[SK]^^V2[SK] 



Kt 



[SK]^-ui[SK] -U2 



(30) 



(31) 



Given 7^, we have R^^ = 0.9R^^ and R^^ = 0.17^. Then 
we can obtain SK^^ and SK^^ from (30) and finally, Kf^ 
and K^^ from (31). This gives the Hill number nn- 

Steady-state approximation using tQSSA 

The tQSSA replaces the free substrate Sp as the slow vari- 
able by the total intact substrate concentration Sp + SpP, 
while retaining the quasi-steady-state assumption for the 
intermediate complexes. Since the extended model has 3 
products {Sp,S and PK), we introduce 2 slow variables 
(total Sp and total PK concentrations). 

Denote Ci = [SpP],C2 = [SpPK],C3 =[SK],C4 = 
[PK] and let X = C2 + C4 to be the total PK concentra- 
tion and A = [Sp] + Ci to be the total Sp concentration, 
respectively. Then it holds that 



P,=[P]+Ci+X, Kt 
St=A+[S]+C2 + C3. 



--[K]+C3+X, 
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The variables A and X satisfy the following equations 
dX 



dt 



dA 



KiiPt - Ci - X){Kt - Cs - X) - K-i{X - C2), 

(32) 

(33) 



— - = -/cid + k2C3 - k3C2^ 

at 

We hypothesize that the intermediate complexes 
Ci, C2, C3 have faster dynamics then the active protein A. 
At steady-state, the intermediate complexes satisfy 

MiCi = {A-Ci){Pt-Ci-X\ (34) 

M3C2 = {A-Ci){X-C2). (35) 

M2C3 = {St-A-C2-C3){Kt-C3-X), (36) 

The simulations of the tQSSA approach are shown in 
Figure 5. 
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